Analytical Models of the Intergalactic Medium and Reionization 
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Reionization is a process whereby hydrogen (and helium) in the Universe is ionized by the radiation 
from first luminous sources. Theoretically, the importance of the reionization lies in its close coupling 
with the formation of first cosmic structures and hence there is considerable effort in modelling 
the process. We give a pedagogic overview of different analytical approaches used for modelling 
reionization. We also discuss different observations related to reionization and show how to use 
them for constraining the reionization history. 



I. INTRODUCTION 

Current models of cosmology ^ indicate that about 
three-fourth of the energy density in our present Universe 
is constituted by "dark energy" which is responsible for 
the current acceleration of the cosmic expansion. The 
next dominant component is the "dark matter" which 
constitutes about 23 per cent of the density. This form 
of matter is coUisionless and interact only gravitationally. 
The baryons constitute only 2 per cent of the total mass. 
The two most abundant elements among the baryons are 
hydrogen and helium. 

Study of reionization mostly concerns with the ioniza- 
tion and thermal history of the baryons (hydrogen and 
helium) in our Universe @j H, ■ Within the frame- 
work of hot Big Bang model, hydrogen formed for the 
first time when the age of the Universe was about 3x10^ 
years, its size being one-thousandth of the present (cor- 
responding to a scale factor a w 0.001 and a redshift 
z = 1/a — 1 ~ 1000). Around this time, the temperature 
of the radiation became low enough w 3 x 10'^ K that 
the photons were not able to ionize the proton-electron 
pair through collisions and hence formation of hydrogen 
(and some amount of helium) could take place Q . The 
epoch at which the protons and the electrons combined 
for the first time to form hydrogen atoms is known as the 
recombination epoch and is well-probed by the Cosmic 
Microwave Background Radiation (CMBR). 

Right after the recombination epoch, the Universe en- 
tered a phase called the "dark ages" where no significant 
radiation sources existed. The hydrogen remained largely 
neutral at this phase. The small inhomogeneities in the 
dark matter density field which were present during the 
recombination epoch started growing via gravitational in- 
stability giving rise to highly nonlinear structures like 
the collapsed haloes 01 • It should, however, be kept in 
mind that most of the baryons at high redshifts do not 
reside within these haloes, they rather reside as diffuse 
gas within the intergalactic space which is known as the 
intergalactic medium (IGM) [1, [§] . 

The collapsed haloes form potential wells whose depth 
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depend on their mass and the baryons (i.e, hydrogen) 
then "fall" in these wells. If the mass of the halo is high 
enough (i.e., the potential well is deep enough), the gas 
will be able to dissipate its energy, cool via atomic or 
molecular transitions and fragment within the halo. This 
produces conditions appropriate for condensation of gas 
and forming stars and galaxies. Once these luminous 
objects form, the era of dark ages can be thought of being 
over. 

The first population of luminous stars and galaxies can 
generate ultraviolet (UV) radiation through nuclear reac- 
tions. In addition to the galaxies, perhaps an early pop- 
ulation of accreting black holes (quasars) also generated 
some amount UV radiation. The UV radiation contains 
photons with energies > 13.6 eV which are then able 
to ionize hydrogen atoms in the surrounding medium, 
a process known as "reionization" . Reionization is thus 
the second major change in the ionization state of hy- 
drogen (and helium) in the Universe (the first being the 
recombination) . 

As per our current understanding [l3, [HI, , reion- 
ization started around the time when first structures 
formed, which is currently believed to be around z « 
20 — 30. In the simplest picture, each source first pro- 
duced an ionized region around it; these regions then 
overlapped and percolated into the IGM. This era is usu- 
ally called the "pre-overlap" phase. The process of over- 
lapping seemed to be completed around z « 6 — 8 at 
which point the neutral hydrogen fraction fell to values 
lower than 10"**. Following that a never-ending "post- 
reionization" (or "post-overlap" ) phase started which im- 
plies that the Universe is largely ionized at present epoch. 
Reionization by UV radiation is also accompanied by 
heating: electron which are released by photoionization 
will deposit an extra energy equivalent to hpV — 13.6 eV 
to the IGM, where v is the frequency of the ionizing pho- 
ton and hp is the Planck constant. This reheating of the 
IGM can expel the gas and/or suppress cooling in the 
low mass haloes - thus, there is a considerable reduction 
in the cosmic star formation right after reionization. In 
addition, the nuclear reactions within the stellar sources 
potentially alter the chemical composition of the medium 
if the star dies via energetic explosion (supernova). This 
can change the nature of star formation at later stages. 

The process of reionization is of immense importance 
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in the study of structure formation since, on one hand, it 
is a direct consequence of the formation of first structures 
and luminous sources while, on the other, it affects sub- 
sequent structure formation. Observationally, the reion- 
ization era represents a phase of the Universe which is yet 
to be probed; the earlier phases {z « 1000) are probed 
by the CMBR while the post-reionization phase {z < 6) 
is probed by various observations based on galaxies, clus- 
ters, quasars and other sources. In addition to the impor- 
tance outlined above, the study of dark ages and cosmic 
reionization has acquired increasing significance over the 
last few years because of the availability of good quality 
data in different areas. 

In this article, we will mainly concentrate on the cur- 
rent status of various analytical and semi-analytical ap- 
proaches which go into modelling reionization. The main 
aim would be to systematically discuss the set of equa- 
tions which are crucial in understanding the process high- 
lighting the major physical processes and assumptions. 
We shall also highlight the relevant observational probes 
at appropriate places. In Section [III we shall give a ped- 
agogic introduction to the basic theoretical formalism for 
studying reionization and IGM in different phases of evo- 
lution. Section |lTT] would be devoted to discussing de- 
tailed modelling of reionization using the formalism de- 
veloped. We shall illustrate on how to constrain the mod- 
els by comparing with a wide variety of available data 
sets. In Section IIVI we shall briefly discuss the current 
numerical simulations and observations related to reion- 
ization. We shall also highlight what to expect in this 
field in near future. 



II. THEORETICAL FORMALISM 

In this section, we discuss the basic theoretical for- 
malism required for modelling reionization of the IGM. 
The main aim here would be to highlight the physical 
processes which are crucial in understanding reionization 
and comparing with observations. In what follows, we 
shall assume that the IGM consists only of hydrogen and 
neglect the presence of helium. It is straightforward to 
include helium into the formalism. 

Essentially, in presence of a ionizing radiation, the evo- 
lution of the mean neutral hydrogen density nni ^148] is 
given by 

nm = -3H{t)nm - Tmnm + Ca{T)nmine (1) 

where overdots denote the total time derivative d/dt, 
H = a/a is the Hubble parameter, Fhi is the photoion- 
ization rate per hydrogen atom, a(T) is the recombina- 
tion rate coefficient and rig represents the mean electron 
density. The first term in the right hand side of equa- 
tion ([1]) corresponds to the dilution in the density be- 
cause of cosmic expansion, the second term corresponds 
to photoionization by the ionizing flux and the third term 
corresponds to recombination of protons and free elec- 
trons into neutral hydrogen. The quantity C is called the 



clumping factor and is defined as 

(nmiX^ie) ("-H>^ 

where the last equality holds for the case when the IGM 
contains only hydrogen (i.e., no helium) and is highly 
ionized, i.e. He ~ tt-hii ~ nn- The clumping factor takes 
into account the fact that the recombination rate in an 
inhomogeneous (clumpy) IGM is higher than a medium 
of uniform density. 

The ionization equation is usually supplemented by the 
evolution of the IGM temperature T, which is given by 

^kin = -2H{t)E^^ + A (3) 

where E'kin = SksTuH is the kinetic energy of the gas 
and A is the net heating rate including all possible heat- 
ing and cooling processes. The first term on the right 
hand side takes into account the adiabatic cooling of the 
gas because of cosmic expansion. 

A. Cosmological radiation transfer 

The equation of radiation transfer, which describes 
propagation of radiation flux through a medium, is writ- 
ten as an evolution equation for the speciflc intensity of 
radiation I^, = I{t, x, v, n) which has dimensions of the 
energy per unit time per unit area per unit solid angle 
per frequency range. It is a function of time and space 
coordinates {t,x), the frequency of radiation i/ and the 
direction of propagation n. The radiation transfer equa- 
tion in a cosmological scenario has the form 

ot a[t) Bu 

c 

= -CK^I„ + 1—^1" (4) 

where is the absorption coefficient and is the emis- 
sivity. The above equation is essentially the Boltzmann 
equation for photons with I^, being directly proportional 
to the phase space distribution function The terms 
on the left hand side of equation (jH) add up to the to- 
tal time derivative of 1,^; in particular, the third term 
corresponds to dilution of the intensity and the fourth 
term accounts for shift of frequency i/ oc a~^ because 
of cosmic expansion. The effect of scattering (which is 
much rarer than absorption in the IGM) can, in princi- 
ple, be included in the Hi, term if required. If the medium 
contains absorbers with number density Uabs each hav- 
ing a cross-section cr^, the absorption coefficient is given 
by Ki, = riabsfj/. The mean free path of photons in the 
medium is given by Ay(t) = K~^{t). 

We define the mean specific intensity by averaging I^, 
over a large volume and over all directions 

f d^x f dn , 
Mt)^ J^ — J -lAt,x,n) (5) 
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Then the spatially and angular-averaged radiation trans- 
fer equation becomes @ 



J. 



dt 



H{t)v 



dv 



47r 



(6) 



where the coefHcients and are now assumed to be 
averaged over the large volume. The quantity is essen- 
tially the energy per unit time per unit area per frequency 
interval per solid angle. 

The integral solution of the above equation along a line 
of sight can be written as [3] 



47r Jo 



r{t.t';v) 



where v' = va{t)/a{t'), v" = va{t) / a{t") and 



T{t,t']v) = c dt"K^,,{t") c 



dt" 



t' K"{t") 



(7) 



(8) 



is the optical depth along the line of sight from t' tot > t' 
Clearly, the intensity at a given epoch is proportional to 
the integrated emissivity with an exponential attenuation 
due to absorption in the medium. The intensity attenu- 
ates by 1/e when the radiation travels a distance equal 
to the mean free path. 

The absorption is "local" when the mean free path of 
photons is much smaller than the horizon size of the Uni- 
verse, i.e., X^{t) <^ c/H{t). In addition, if we also assume 
that the emissivity does not evolve significantly over 
the small time interval A/c, then the specific intensity is 
related to the emissivity through a simple form [l^ [l6j 



47r 



(9) 



Note that in the case of local absorption, J, HJ <C ce. 
In this approximation, the background intensity depends 
only on the instantaneous value of the emissivity (and 
not its history) because all the photons are absorbed 
shortly after being emitted (unless the sources evolve 
synchronously over a timescale much shorter than the 
Hubble time). We shall discuss later in Section Hi BI that 
this is a useful approximation for the IGM for redshifts 
^>3. 



B. Post-reionization epoch 

Let is first study the radiation transfer in the post- 
reionization epoch. Compared to the pre-overlap era, this 
epoch is much easier to study because the IGM can be 
treated as a highly ionized single-phase medium (whereas 
during the pre-overlap era, one is looking into two dis- 
tinct phases - ionized and neutral). The optical depth 
can be written as 



t(z, z'; v) = c 



dz" 



(1 + z")H{z") 



nm{z")a^Uv") (10) 



where dabs (i^) is the total absorption cross section of neu- 
tral hydrogen and we have changed the time coordinate 
to the redshift z. Various processes can, in principle, 
contribute to crabs(i^)7 most dominant being the resonant 
Lyman series absorption corresponding to excitation of 
hydrogen atoms from the ground state to higher ones (Is 
— > np) and the continuum absorption of photons above 
the ionization threshold via photoionization process. Let 
us treat each of them separately in the following: 



1. Resonant Lyman series absorption 

The Lyman series absorption arises from the electronic 
excitation of neutral hydrogen atoms from the Is ground 
state to higher ones. The most dominant of these are the 
Lya (Is 2p, rest wavelength Aq, « 1216A ) and Ly/3 (Is 

3p, A/3 ~ 1206^1 ) transitions, and hence they are the 
most relevant ones as far as observations are concerned. 
For simplicity, we shall present results for the Lya ab- 
sorption only, the others can be calculated in identical 
manner. The Lya absorption cross section is given by 



0-abs(,l^) 



1 



(11) 



where Va — c/Aq is the resonant frequency of transition, 
a a = 4.45 X 10~^^cm^ is the cross section at I'aipha and 
y is a function which determines the profile of the ab- 
sorption line. It is called the Voigt profile function and 
is a convolution of the Lorentzian shape for the natu- 
ral broadening and the Gaussian shape for the thermal 
broadening. For the purpose of this article, it is sufh- 
cient to note that y is a sharply peaked function about 
v/va = 1; for most our discussion, we shall take it to be 
a Dirac-delta function V{i'li'a ~ 1) = ^Diyjva — 1). 

The optical depth between the redshifts z and z' is 
then given by 



t(z, z'; v) = t(zq) 



■H(z„) 



1 + Za 



nYa{za)\ 
-(1 + -) 



(12) 



If we put this into equation ([7]), we see that the Lya 
absorption at a redshift Zq reduces the specific intensity 
observed at z at a frequency v aiX + z) j (\ ^ Za) by a fac- 
tor e"'^^^"^. The value of T(za) along a given line of sight 
would depend upon the distribution of nYa{za). How- 
ever, we would mostly be interested in the mean value 
of specific intensity averaged over a number of lines of 
sight. The corresponding reduction can be described by 
a line-of-sight-averaged optical depth 



cff(^c.) ^ 



LOS 



(13) 



where ( )los denotes averaging over lines of sight. The 
quantity Toff is usually known as the "effective optical 
depth" . 
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FIG. 1: The eflective optical depth of Lya absorption as func- 
tion of redshift z. The points with errorbars represent the 
observational data. The dashed curves, from top to bottom, 
represent the predictions for a uniform IGM with neutral hy- 
drogen fraction xhi ~ (3, 1,0.3) x 10"'', respectively. 



Theoretically, the value of Teg can be calculated if 
we know the distribution of optical depth P(t) [which 
can be calculated from the neutral hydrogen distribution 
P{nm)]- 



-Toff (z) 



dr P{t; z) e 



(14) 



Of course, one requires detailed understanding of the evo- 
lution of the baryonic density field to model the distri- 
bution P(nHi). This has been addressed by a series of 
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23 241 and 
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analytical studies [l7l ll8l ll9l l20l |21| 
numerical simulations [2^ [2^, 123, [2^, |28 
which we shall avoid discussing here. However, we can 
still make some inference assuming the distribution is 
uniform, i.e., Toff(z) = caanmiz)/ H{z). If we define 
the neutral hydrogen fraction to be xhi = nm/nn, then 
we can calculate oc a;Hi given a set of cosmological 
parameters (which would uniquely determine H{z) and 
uh). 

Observationally, Tcs can be determined by looking at 
the spectra of bright sources like quasars at high red- 
shifts. These spectra show a series of absorption fea- 
tures at frequencies larger than the Lya frequency in the 
quasar rest frame. Since one has a good knowledge of 
the unabsorbed quasar spectra (from looking at nearby 
quasars and also having some understanding about the 
physical processes), one can calculate the amount of ab- 
sorption happening because of the intervening IGM be- 
tween the quasar and the observer; this absorption, av- 
eraged over numerous lines of sight, is essentially the 
quantity e"'^"^'. This has been done to quite high red- 
shifts z ~ 6.5 and the values of Tes observed are shown 
as points with errorbars [s^, UH] in Figure [TJ 

To understand what these values imply, we have plot- 
ted with dashed lines the calculated value of for a 
uniform IGM assuming three values of xhi = (3, 1, 0.3) x 



10~^ from top to bottom. This immediately tells us that 
the fraction of neutral hydrogen has to be ^ 10~^ in 
order to reproduce the observed values of Tcs- In fact, 
had xhi been slightly (say ~ 10~^) higher, one would 
have obtained Tcff much higher than unity (~ 10 — 100) 
and hence the flux from the quasar would be completely 
absorbed. If that were the case, it would show up as a 
absorption "trough" at frequencies larger than the rest 
frame Lya frequency. In reality, a considerable amount 
of transmitted flux is found at these frequencies along- 
with a series of absorption features arising from the Lya 
transition of residual neutral hydrogen. These absorp- 
tion signatures are known as the "Lya" forest and are 
powerful probes of the neutral hydrogen distribution in 
the IGM at z < 6 m. 

The absence of a absorption trough is a direct proof 
of the fact that hydrogen is completely reionized in the 
diffuse IGM at redshifts z ^ 6. This is known as the 
Gunn-Peterson effect ^33]. Note that the actual inferred 
value of Xm might be slightly different if one models with 
an appropriate density distribution, however, the basic 
conclusion remains unchanged. 

We can also see from Figure [T] that for quasars at red- 
shifts z ^ 6, the observed value of TeS ^ 5; this would 
imply an attenuation > 0.01 and hence one actually ob- 
serves absorption troughs as predicted by Gunn-Peterson 
effect l35i, ^] . Unfortunately, finding such troughs does 
not necessarily imply that the IGM is highly neutral as 
even a xhi ^ lO""* could be sufficient to absorb all the 
flux. However, one can use much detailed modelling to 
improve the constraint, which we shall discuss later in 
Section UnBl 

The values of Tcff < 1 at z < 4 means that the diffuse 
IGM is highly transparent (also called optically thin) to 
Lya photons. Only about ~ 10% of the Lya photons 
are absorbed, mostly within the high density regions. 
These high density systems are often modelled as a set 
of discrete absorbers of some size. If we consider an ab- 
sorber having neutral hydrogen density nni and a size 
L ^ c/H{z) along the line of sight at a redshift Zabs, the 
optical depth is given by 



t(z, z'] v) = A^HI a'abs(i^abs) 



(15) 



where i^abs = ^^(1 + ^abs)/(l + z) and TVhi = nni L is the 
column density of neutral hydrogen within the absorber. 
Hence, each absorber reduces the specific intesnity by a 
factor e~^Hio-abi,(!^abB)^ If we assume that the absorbers are 
Poisson-distributed, then it is straightforward to show 
that the effective optical depth is given by ^] 



roff(z, z'; u) 
dz 



°° rfiN 



dz"dNm 



[1-e 



](16) 



where [d^N / dzdN^a] dz diVni is the number of absorbers 
within (z, z -|- dz) having column densities in the range 
{Nm,Nm + dNm)- 



5 



In case of Lya resonant absorption, we can use the 
cross section in pT|) to calculate Teff. Since F is a func- 
tion which is sharply peaked around vjva — 1, we can 
approxiamte the above integral as 



T,s{z, z'- u) = / diVHi ."^.^^ W^{N^i) 

K Jo dzaONm 



where 

Wa^iNm) = / dA" 



I _ g-iVHlcrcV(A„/A"-l) 



(17) 



(18) 



is called the "equivalent width" of the absorber. 



2. Continuum absorption 

In case of continuum absorption of radiation by pho- 
toionization, the cross section is given by 



o'abs(i^) = crm{v)Q{i^ - mi) 



(19) 



where (Thi(i^) is the photoionization cross section and Q 
is the Heaviside step function taking into account that 
only photons with frequencies v > vm = 13.6 eV /hp 
would be absorbed by the photoionization process. The 
exact form of (Jbi{v) is rather complicated, however one 
can approximate it by a power-law of the form UYniv) = 
<yo{v/vm)~^ where ctq = 6.3 x 10~^*cm^. 

Since ctq ~ ctq, one can show that the absorption due 
to a diffuse IGM in this case too is negligibly small. The 
only significant absorption can be seen in very high den- 
sity regions which have a large fraction of their hydrogen 
in neutral form. In that case, we can use the relations 
obtained for a set of Poisson-distributed absorbers in the 
vase of resonant transition. We essentially have an opti- 
cal depth of the form (|17p . and the corresponding mean 
free path of ionizing photons due to these discrete ab- 
sorbers is found to be [13, SH 



H{z){l + z) 



diV, 



HI 



92 AT 
dzdN] 



_ Nuirr-aiiv) 



HI 



(20) 



At this point, let us introduce the concept of Lyman- 
limit systems which have column densities -/Vhi > cr^^ — 
1.6 X 10^^cm~2; these absorbers contribute a optical 
depth of unity to the ionizing photons. The average dis- 
tance between these systems is given by 



A 



LLS 



Hiz)il 



dAT] 



LLS 



dz 



(21) 



where dA'LLs/d^ is the redshift distribution of the 
Lyman-limit systems 



d^LLS 

dz 



dTVi 



HI 



l/o-HlCl'HI 



dzdNi 



HI 



(22) 



For the observed distribution d^N/dzdNm cx iVni 
one can show from equations (PO)) and (PT|) that the mean 
free path is related to the distance and redshift distribu- 
tion Lyman-limit systems as 

V^i7(z)(l + z) diVLLs/dz- ^^-"^ 

The redshift distribution of Lyman-limit systems 
dA^LLg/dz is a quantity which has been measured for 
2 < z < 4.5 by observations of quasar absopriton spectra. 
Though the observational constraints are poor, one can 
stih obtainavalue dA^LLs/dz « 0.3{l+z)-^-^^ which 
in turn gives the mean free path as X„^^/[c/ H{z)] w 
0.1[(1 -I- z)/4]~2-^^ Hence the mean free path of ionizing 
photons is much smaller than the horizon size for z > 3, 
which implies that we can use the local absorption ap- 
proximation at these redshifts. 

We can summarise the main results of this section as: 
the post-reionization epoch is characterized by a highly 
ionized IGM as observed by the quasar absorption spec- 
tra. The IGM is largely transparent to ionizing photons 
at these redshifts. However, there exist regions with high 
column densities (A'hi > lO^^cm"^) which are optically 
thick to the ionizing radiation; these regions determine 
the photon mean free path. We shall see later how to 
use this information to obatin an improved model of the 
IGM. 



C. Pre-overlap epoch 

We now turn our attention towards the IGM in the 
pre-overlap era. In this era, the overlap of individual 
ionized regions is not complete and hence the IGM is 
partially ionized. So the radiation transfer equation has 
to be modified to account for the multi-phase nature of 
the IGM. 

Let us define the volume filling factor of ionized regions 
to be Qhii; this the fraction of volume that is ionized and 
reionization is said to be complete when Qhii — 1- Next, 
note that the number density of photons present in the 
background flux is 



nj{t) 



47r 
c 



Ju_ 



(24) 



Since there is no ionizing flux within the neutral regions 
(otherwise they would not remain neutral) , the photoion- 
ization rate per hydrogen atom within the ionized (HII) 
regions is 



HI 



1 f°° J,, 

- — 47r / dv—^GYa{v) 

/hII Jy-ai n-pV 



(25) 



where the factor Qjjjj accounts for the fact that the ra- 
diation is limited to a fraction of the total volume. The 
emmission rate of ionizing photons per unit volume from 
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sources of emissivity is 



D. Reionization of the inhomogenous IGM 



(26) 



Then the equation of radiation transfer then becomes 

[mil 



-3H{t)nj - H{t) 



47r .L 



riph 



dQ 



HII 



HII" 



dt 



(27) 



where n^j and the number densities of neutral 

and ionized hydrogen within the HII regions, respectively. 
The first term in the right hand side of equation P7|) 
corresponds to the dilution in density due to cosmic ex- 
pansion while the second term accounts for the loss of 
ionizing radiation because of a photon being redshifted 
below the ionization edge of hydrogen z^hi- The third 
term is essentially the source of ionizing photons. The 
fourth term accounts for the loss of photons in ionizing 
the residual neutral hydrogen within the ionized regions. 
The fifth term, which is only relevant for the pre-overlap 
stages, accounts for the photons which ionize hydrogen 
for the first time and hence increase the filling factor 
Qhii- For Qhii — 1, equation (P7|) reduces to that for 
the post-reionization phase. 

If we now assume that the photons are absorbed lo- 
cally, then J, HJ ce and we is essentially given by 
equation We can then ignore terms containing J and 
nj in equation (|27p . This gives a equation describing the 
evolution of the filling factor Qhii 



dQ 



HII 



dt 



— Qhii! hi 

"hii 



(28) 



If we further assume photoionization equilibrium within 
the ionized region d{n^^a^) / dt — > 0, then we have from 
equation ^ "hiFhi — Ca(T)n\^jjn}} and the evolution 
of Qhii can be written in the form [1^ 



dQ 



HII 



dt 



'HII 



Qhii Ca{T)n]} 



(29) 



In this description, reionization is complete when Qhii — 
1 and equation ([29]) cannot be evolved further on. Clearly 
the assumptions of local absorption and photoionization 
equilibrium (both of which are reasonably accurate) has 
given us a equation which can be solved once we have 
a model for estimating and C. Of course, there is a 
dependence of the recombination rate coefficient a{T) on 
temperature, however that dependence is often ignored 
while studying the volume filling factor. In case one is in- 
terested in temperature evolution, one has to solve equa- 
tion ([3]) taking into account all the heating and cooling 
processes in the IGM, in particular, the photoheating by 
ionizing photons whose rate is given by 



r 



ph,HI 



An 



di^ 



hpV 



(30) 



The description of reionization in the previous section 
is not adequate as it does not take into account the in- 
homogeneities in the IGM appropriately (except for a 
clumping factor C in the effective recombination rate). 
To see this, consider the post-reionization phase where we 
know from observations that there exists regions of high 
density which are neutral; these regions are being grad- 
ually ionized and hence one would ideally like to write a 
equation similar to (|29p for studying the post-reionization 
phase. Since the ionization state depends on the density, 
one should have to account for the density distribution 
of the IGM. 

In order to proceed, first note that the volume fill- 
ing factor may not be the appropriate quantity to study 
for evolution of reionization because most of the pho- 
tons are consumed in regions with high densities (which 
might be occupying a small fraction of volume). In other 
words, if we neglect recomibation for the moment, we 
have from equation (|^^ that the volume filling factor 
Qhii = / hph/nn — npi^/nn] however, in reality the 
photon to hydrogen ratio should be equal to the ionized 
mass fraction i^nfj, i.e., Upi^/njj — i^HH- Hence, we must 
replace the volume filling factor by the mass filling factor 
in the description of the previous section, in particular 
equation P5|) should have the form 



dt 



'HII 



(31) 



One can relate Fj^fj to the IGM density distribution by 
using the fact that regions of lower densities will be ion- 
ized first, and high-density regions will remain neutral for 
a longer time. The main reason for this is that the recom- 
bination rate (which is oc n'j^) is higher in high-density 
regions where dense gas becomes neutral very quickly. If 
we assume that hydrogen in all regions with overdensities 
A < Ahh is ionized while the rest is neutral, then the 
mass ionized fraction is clearly (l6| 



F-i^^^^''^(Ahh) 



dA A P(A) (32) 



where -P(A) is the (volume-weighted) density distribu- 
tion of the IGM. The term describing the effective re- 
combination rate gets constribution only from the low 
desnity regions (high density neutral regions do not con- 
tribute) and is then given by 



At 



a{T)nl' / dA A^ P(A) = a{T)nl' R{A^u) (33) 
Jo 

The evolution for the mass ionized fraction is then 

dF-'^(AHii) hph , , , H /g^\ 
= " R{Ami)a{T)n^ (34) 



'HII 



The evolution equation essentially tracks the evolution of 
Ahh which rises as i^^^(AHii) increases with time (i.e.. 
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more and more high density regions are getting ionized). 
Since the mean free path is determined by the high den- 
sity regions, one should be able to relate it to the value of 
Ahii 16]. It is clear that a photon will be able to travel 
through the low density ionized volume 



Fy(AHii) = / ™dA P(A) 



(35) 



before being absorbed. When a very high fraction of 
volume is ionized, one can assume that the fraction of 
volume filled up by the high density regions is 1 — Fy, 
hence their size is proportional to (1 — FvY^^, and the 
separation between them along a random line of sight 
will be proportional to (1 — Fv)~^^^, which, in turn, will 
determine the mean free path. Then one has 



A^(a) = 



An 



[1-Fv(Ahii)]2/3 



(36) 



where we can fix Ao by comparing with low rcdshift ob- 
servations like the distribution of Lyman-limit systems 
[equation ([23]) ]. 

The situation is slightly more complicated when the 
ionized regions are in the pre-overlap stage. At this stage, 
a volume fraction 1 — Qhii of the universe is completely 
neutral (irrespective of the density), while the remain- 
ing Qhii fraction of the volume is occupied by ionized 
regions. However, within this ionized volume, the high 
density regions (with A > Ahii) will still be neutral. 
Once Qhii becomes unity, all regions with A < Ahii are 
ionized and the rest are neutral; this can be thought of as 
the end of reionization. The generalization of equation 
P4p . appropriate for this description is given by 



equation 



d[QHii-FM(AHii)] _ np]^{z) 



dt 



QmiaR{T)neR{Ami) 



mil 



(37) 

Note that there are two unknowns Qhii and F]\i{A}iii) 
in equation (|37p which is impossible to solve it with- 
out more assumptions. One assumption which is usu- 
ally made is that Ahii does not evolve significantly with 
time in the pre-overlap stage, i.e., it is equal to a crit- 
ical value Ac. This critical density is determined from 
the the mean separation of the ionizing sources. To have 
some idea about the value of Ac, two arguments have 
been put forward in the literature: In the first, it is ar- 
gued that Ac is determined by the distribution of sources 
[lOlIT^. When the sources are very numerous, every low- 
density region (void) can be ionized by sources located 
at the edges, and hence the overlap of ionized regions 
can occur (i.e., Qhii approaches unity) when Ac ^ 1 is 
the characteristic overdensity of the thin walls separat- 
ing the voids. For rare and luminous sources, the mean 
separation is much larger and hence the value of Ac has 
to be higher before Qhii can be close to unity. In the 
second approach, it is assumed that the mean free path 
is determined by the distance between collapsed objects 
(which manifest themselves as Lyman-limit systems) and 



hence Ac should be similar to the typical overdensities 
near the boundaries of the collapsed haloes [i^ [4^ . It 
usually turns out to be ^ 50 — 60 depending on the den- 
sity profile of the halo. Interestingly, results do not vary 
considerably as Ac is varied from ~ 10 to ^ 100. Once 
Ac is fixed, one can follow the evolution of Qhii until it 
becomes unity. Following that, we enter the post-overlap 
stage, where the situation is well-described by equation 

Of course, the above descrition is also not fully ade- 
quate as there will be a dependence on how far the high 
density region is from an ionizing source. A dense region 
which is very close to an ionizing source will be ionized 
quite early compared to, say, a low-density region which 
is far away from luminous sources. However, it has been 
found that the above description gives a reasonable an- 
alytical description of the reionization process, particu- 
larly for the post-reionization phase. The main advan- 
tages in this approach are (i) it takes into account both 
the pre-overlap and post-overlap phases under a single 
formalism, (ii) once we have some form for the IGM den- 
sity distribution P(A), we can calculate the clumping fac- 
tor and the effective recombination rate self-consistently 
without introducing any extra parameter; in addition we 
can also compute the mean free path using one single 
parameter (Aq, which can be fixed by comparing with 
low-redshift observations). 



III. MODELLING OF REIONIZATION 

Given the formalism we have outlined in the previous 
section, we can now go forward and discuss some other 
details involved in modelling reionization. 

A. Reionization sources 

The main uncertainty in any reionization model is to 
identify the sources. The most natural sources which 
have been observed to produce ionizing photons are the 
star-forming galaxies and quasars. Among these, the 
quasar population is seen to decrease rapidly at z > 3 
and there is still no evidence of a significant population 
at higher redshifts. Hence, the most common sources 
studied in this area are the galaxies. 

The subject area of formation of galaxies is quite in- 
volved in itself dealing with formation of non-linear struc- 
tures (haloes and filaments), gas cooling and generation 
of radiation from stars. 



1. Mass function of collapsed haloes 

The crucial ingredient for galaxy formation is the col- 
lapse and virialization of dark matter haloes. This can 
be adequately described by the Press-Schechter formal- 
ism for most purposes. It can be shown that the number 
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density of collapsed objects per unit comoving volume 
(which is physical volume divided by a^) within a mass 
range (A/, M + dM) at an epoch t is given by [i^l 



dn{M, t) 
dM 



dM 



2 ./,p^dlna(M)^ (38) 
M dlnM M ' ^ ' 



where Pm is the comoving density of dark matter, (7(M) 
is defined as the rms mass fluctuation at a mass scale M 
at z = 0, 1/ = Sc/[D{t)aiM)], D{t) is the growth factor 
for linear dark matter perturbations and 5c is the critical 
overdensity for collapse, usually taken to be equal to 1.69 
for a matter-dominated flat universe = !)• This for- 
malism can be extended to calculate the comoving num- 
ber density of collapsed objects having mass in the range 
(M, M -|- dM), which are formed within the time inter- 
val (ifgrm ) iform + diform) ^ud obscrvcd at a later time t 
is given by 



-aMdtform = 



dMdtto 



dMdtio 

^ Psuiv 



for 



dMdiform 



(39) 



where d'^n{M, iform)/9M9iform|form is the formation rate 
of haloes at tform and Psurv[i Inform] is their survival prob- 
ability till time t. Assuming that haloes are destroyed 
only when they merge to a halo of higher mass, both 
these quantities can be calculated from the merger rates 
of haloes. The merger rates can be calculated using de- 
tailed properties of gaussian random field. The quantities 
can also be calculated in a more simplistic manner by as- 
suming that the merger probability is scale invariant; in 
that case (H 



d^n{M,t) 



dMdt 



d^n{M,t) dn{M,t)D 



form 



dMdt 

dn{M,t) 2 D 
dM ^ D 



dM D 



and 



tltfo 



^surv l'^ t-jorm 



D{t form) 



(40) 



(41) 



2. Star formation rate 

If these dark matter haloes are massive enough to form 
huge potential wells, the baryonic gas will simply fall 
into those wells. As the gas begins to settle into the 
dark matter haloes, mergers will heat it up to the virial 
temperature via shocks. However, to form galaxies, the 
gas has to dissipate its thermal energy and cool. If the 
gas contains only atomic hydrogen, it is unable to cool 
at temperatures lower than 10* K because the atomic 
hydrogen recombines and cannot be ionized by collisions. 
The gas can cool effectively for much lower temperatures 



in the presence of molecules - however, it not at all clear 
whether there are sufficient amount of molecules present 
in the gas at high redshifts [i^, [s^, [EH . Hence the lower 
mass cutoff for the haloes which can host star formation 
will be decided by the cooling efficiency of the baryons. 

Let M^,(M,t,t[oim) denote the rate of star formation 
at time t within a halo of mass M which has formed at 
^form- Then we can write the cosmic SFR per unit volume 
at a time i. 



P*it) = I dtform I dM'M,{M\t,tfo 

d^n{M' ,t]tform) 
dM'dtlorm 



(42) 



where the is included to covert from comoving to 
physical volume. The lower mass cutoff Mmin{t) at a 
given epoch is decided by the cooling criteria as explained 
above. However, once reionization starts and regions are 
reheated by photoheating, the value of Mmin{t) is set 
by the photoionization temperature ~ 10** K. This can 
further suppress star formation in low mass haloes and is 
known as radiative feedback. We shall discuss this later 
in this section. 

The form of M*(M, iform) contains information 
about various cooling and star-forming processes and 
hence is quite complex to deal with. It can be obtained 
from semi-analytical modelling of galaxy formation (s^ 
or constrained from observations of galaxy luminosity 
function [53| . A very simple assumption that is usually 
made for modelling reionization is that the duration of 
star formation is much less than the Hubble time H~'^{t) 
which is motivated by the fact that most of the ionizing 
radiation is produced by hot stars which have shorter life- 
time. In that case, M*(M, i, tform) can be approximated 
as 

M,{M,t,tform) ~ MJoit-tform) = /* ^Mfo (i- ^form) 

Pm 

(43) 

where pb/pmM is the mass of baryons within the halo 
and /* is the fraction of baryonic mass which has been 
converted into stars. The cosmic SFR per comoving vol- 
ume is then 



P*(i) 



1 Pb 
a^{t) Pm 



dM' /, Af 



, d'^n{M',t) 



M„,„(t) 



(44) 

One should keep in mind that many details of the star- 
formation process has been encoded within a single pa- 
rameter This should, in principle, be a function of 
both halo mass M and time t. However, it is not clear at 
all what the exact dependencies should be. Given such 
uncertainties, it is usual to take it as a constant. 

In addition, one finds that the merger rate of haloes 
at high redshifts is much less than the formation rate 
(which follows if 1/ > 1), hence d^n{M,t)/dMdt\form ~ 
d'^n{M, t)/dMdt. Then, one can write the SFR in terms 
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of the fraction of collapsed mass in haloes more massive 
than M„iin(t) 



1 

fcoiiit) - — / dM 

Pm JAf„i„(t) 



, ^^^, dn{M',t) 



erfc 



as 



P*{t) = /. 



V2 D{t)aiM^in) 
Pb d/coii(t) 



(45) 



(46) 



where ni, is the total baryonic number density in the IGM 
(equal to riH if we ignore the presence of helium) and 



-^ion — f^fcsc^p I 



dAL 



(49) 



is the number of photons entering the IGM per baryon 
in collapsed objects [6l|. In case there are more than 
one population of stars, one has to use different values of 
iVion for the different populations. 



4- Feedback processes 



3. Production of ionizing photons 

Given the SFR, we can calculate the emissivity of 
galaxies, or equivalently the rate of ionizing photons in 
the IGM per unit volume per unit frequency range: 



it) = u. 



dN, 
dNL 



P*{t) 



(47) 



where dN^/dM^ gives the number of photons emitted 
per frequency range per unit mass of stars and /esc is the 
fraction of ionizing photons which escape from the star 
forming haloes into the IGM. The emissivity is simply 

Given the spectra of stars of different masses in a 
galaxy, and their Initial Mass Function (IMF), dNi^/dM^, 
can be computed in a straightforward way using "pop- 
ulation synthesis" codes 54 , [s^ . The IMF and spectra 



will depend on the details of star formation (burst forma- 
tion or continuous) and metallicity. In fact, it is possible 
that there are more than one population of stellar sources 
which have different dN^/dM^,. For example, there are 
strong indications, both from numerical simulations and 
analytical arguments [56j . that the first generation stars 
were metal-free, and hence massive, with a very different 
kind of IMF and spectra than the stars we observe today 
[13]; they are known as the Population III (PopIII) stars. 

A fraction of photons produced in a galaxy would 
be consumed in ionizing the neutral matter within the 
galaxy itself. Hence only a fraction of photons escape 
into the IGM and is available for reionization, which is 
encoded in the parameter ffisn- This parameter is again 
not very well modelled [5^, [53. [60. '6r| and its observed 
value is also quite uncertain [6^,163, 64, 65]; typical values 
assumed are ~ 0.1. It is most likely /osc, like is also 
a function of halo mass and the time of halo formation, 
however since the dependences are not well understood, 
it is taken to be a constant. 

The total number of ionizing photons is then obtained 
by integrating the above quantity over all energies above 
the ionization threshold: 



it) 



d/coii(t) 
dt 



(48) 



The moment there is formation of stars and other lu- 
minous bodies, they start to affect the subsequent forma- 
tion of structures - this is known as feedback [g^l- The 
process is intrinsically non-linear and hence quite com- 
plex to model. The feedback processes can be categorized 
roughly into three categories. 

The first of them is the radiative feedback which is 
associated with the radiation from first stars which can 
heat up the medium and can photoionize atoms and/or 
photodissociate molecules. Once the first galaxies form 
stars, their radiation will ionize and heat the surrounding 
medium, increasing the mass scale (often referred to as 
the filtering mass [22| ) above which baryons can collapse 
into haloes within those regions. The minimum mass 
of haloes which are able to cool is thus much higher in 
ionized regions than in neutral ones. Since the IGM is 
multi-phase in the pre-overlap phase, one needs to take 
into account the heating of the ionized regions right from 
the beginning of reionization. In principle, this can be 
done self-consistently from the evolution of the tempera- 
ture of the ionized region in equation ([3]). 

The low mass haloes can be subjected to mechanical 
feedback too, which is mostly due to energy injection via 
supernova explosion and winds. This can expel the gas 
from the halo and suppress star formation. As in the case 
for radiative feedback, one can parametrize this through 
the minimum mass parameter M^i^{t). 

Finally, we also have chemical feedback where the stars 
expel metals into the medium and hence change the 
chemical composition. This would mean that the sub- 
sequent formation of stars could be in a completely dif- 
ferent environment and hence the nature of stars would 
be highly different. 



5. Quasars 

Besides the stellar sources, a population of accreting 
black holes or quasars are also known to produce signif- 
icant amount of ionizing radiation. Hence it is possible 
that they too have contributed to reionization. The frac- 
tion of their contribution would depend on the number 
of quasars produced at a particular redshift. Observa- 
tionally, the luminosity function of quasars is quite well- 
probed till a redshift z ~ 6 It turns out that the 
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population peaks around z « 3 and decreases for higher 
redshifts. Hence their contribution at higher redshifts is 
highly debated. 

The difference between reionization by stellar sources 
and quasars lie in the fact that quasars produce signifi- 
cant number of high energy photons compared to stars. 
This would imply that quasars can contribute signifi- 
cantly to double-reionization of helium (which requires 
photons with energies > 54.4 eV, not seen in galaxies). 
In addition, quasars produce significant amount of X-ray 
radiation. Since the absorption cross section of neutral 
hydrogen varies with frequency approximately as h'~^, 
the mean free path for photons with high energies would 
be very large. A simple calculation will show that for 
photons with energies above 100-200 eV, the mean free 
path would be larger than the typical separation between 
collapsed structures [1^ (the details would depend upon 
the redshift and exact description of collapsed haloes). 
These photons would not be associated with any partic- 
ular source at the moment when they are absorbed, and 
thus would ionize the IGM in a more homogeneous man- 
ner (as opposed to the overlapping bubble picture for UV 
sources) . 



lution of Lyman- limit absorption systems dA^LLS / d-z [zH , 
(ii) the effective optical depths Tes = — In F for Lya and 
Ly/3 absorption in the IGM [s^l, (iii) electron scattering 
optical depth t^i — utc J dt rie (where ax is the Thomson 
scattering cross section) as measured from CMBR experi- 
ments [72| , (iv) temperature of the mean intergalactic gas 
[t^, (v) cosmic star formation history 74], and (vi) 
source number counts at z « 10 from NICMOS HUDF 

The data constrain the reionization scenario quite 
tightly. We find that hydrogen reionization starts at 
z ~ 15 driven by metal-free (PopIII) stars, and it is 90 
per cent complete by z « 8. After a rapid initial phase, 
the growth of the volume filled by ionized regions slows 
down at z ^ 10 due to the combined action of chemical 
and radiative feedback, making reionization a consider- 
ably extended process completing only at z « 6. The 
number of photons per hydrogen at the end of reioniza- 
tion at z w 6 is only a few, which implies that reionization 
occurred in a "photon-starved" manner [76| . 



C. Clustering of sources 



B. Illustration of a semi-analytical model 

The physics described above in the preceding sections 
can all be combined to construct a semi-analytical model 
for studying the thermal and ionization history of the 
IGM. We shall give an explicit example of one such model 
[T^ . [i^ whose main features are: The model accounts 
for IGM inhomogeneities by adopting a lognormal dis- 
tribution for P{A); reionization is said to be complete 
once all the low-density regions (say, with overdensities 
A < Ac ~ 60) are ionized. The ionization and thermal 
histories of neutral, HII and Helll regions are followed 
simultaneously and self-consistently, treating the IGM 
as a multi-phase medium. Three types of reionization 
sources have been assumed: (i) metal-free (i.e. PopIII) 
stars having a Salpeter IMF in the mass range 1 — lOOAf©: 
they dominate the photoionization rate at high redshifts; 
(ii) PopII stars with sub-solar metallicities also having a 
Salpeter IMF in the mass range 1 — IOOMq; (iii) quasars, 
which are significant sources of hard photons at z < 6; 
they have negligible effects on the IGM at higher red- 
shifts. 

As discussed earlier, reionization is accompanied by 
various feedback processes, which can affect subsequent 
star formation. In this model, radiative feedback is com- 
puted self-consistently from the evolution of the ther- 
mal properties of the IGM. Furthermore, the chemical 
feedback inducing the PopIII — > PopII transition is im- 
plemented using a merger-tree "genetic" approach which 
determines the termination of PopIII star formation in a 
metal-enriched halo [70l |. 

The predictions of the model are compared with a wide 
range of observational data sets, namely, (i) redshift evo- 



The formalism described till now works quite well for 
studying global properties of reionization. However, it is 
not adequate for studying the details of overlap of ion- 
ized regions in the pre-reionization epoch. The shapes 
and distribution of ionized regions are crucially depen- 
dent on the distribution of sources, in particular, their 
clustering and also on the density structure of the IGM. 
For example, if the sources are highly clustered, then it 
is expected that the overlap of the ionized bubble for 
nearby sources would be much earlier than what is ex- 
pected from a random distribution of sources. 

Modelling of reionization to such details have become 
important because of the 21cm observations which are 
likely to be available in near future. Before discussing the 
theoretical model, let us briefly outline the motivation 
behind 21cm experiments. 



1. 21 cm observations 

Perhaps the most promising prospect of detecting the 
fluctuations in the neutral hydrogen density during the 
(pre-)reionization era is through the 21 cm emission ex- 
periment s ITTI . some of which are already taking data 
(GMRT dU 21CMA [T50l), and some are expected in 
future (MWA [l5l|, LOFAR 152J, SKA [153.] )■ The ba- 
sic principle which is central to these experiments is the 
neutral hydrogen hyperfine transition line at a rest wave- 
length of 21 cm. This line, when redshifted, is observable 
in radio frequencies (^ 150 MHz for z ^ 10) as a bright- 
ness temperature: 

dTb(z,n)^ — — „ (50) 



l-fz 
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FIG. 2: Comparison of analytical model predictions with observations for the best-fit model. The different panels indicate: (a) 
The volume-averaged neutral hydrogen fraction xm, with observational lower limit from quasar absorption lines at z = 6 and 
upper limit from Lya emitters at z = 6.5 (shown with arrows). In addition, the ionized fraction Xe is shown by the dashed line, 
(b) SFR p, for difi'erent stellar populations, (c) The number of source counts above a given redshift, with the observational 
upper limit from NICMOS HUDF shown by the arrow. The contribution to the source count is zero at low redshifts because 
of the J-dropout selection criterion, (d) Electron scattering optical depth, with observational constraint from WMAP 3-year 
data release, (e) Lya effective optical depth, (f) Ly/3 effective optical depth, (g) Evolution of Lyman-limit systems, (h) 
Photoionization rates for hydrogen, with estimates from numerical simulations (shown by points with error-bars; jSlJ). (i) 
Temperature of the mean density IGM. 



where Ts is the spin temperature of the gas, Tcmb = 
2.76(1 -I- z) K is the CMBR temperature, Aio is the Ein- 
stein coefficient and vq — 1420 MHz is the rest frequency 
of the hyperfine Une. 

The observabihty of this brightness temperature 
against the CMBR background will depend on the rel- 
ative values of T5 and Tcmb- Depending on which pro- 
cesses dominate at different epochs, T5 will couple ei- 
ther to radiation (Tcmb) or to matter (determined by 
the kinetic temperature Tk) Almost in all models 

of reionization, the most interesting phase for observing 
the 21 cm radiation is 6 ^ z ^ 20. This is the phase 
where the IGM is suitably heated to temperatures much 
higher than CMBR (mostly due to X-ray heating (TOj l 
thus making it observable in emission. In that case, we 
have 5Tb oc nm/riH, which means that the observations 
would directly probe the neutral hydrogen density in the 
Universe. Furthermore, this is the era when the bubble- 
overlapping phase is most active, and there is substantial 
neutral hydrogen to generate a strong enough signal. At 
low redshifts, after the IGM is reionized, tt-hi falls by or- 
ders of magnitude and the 21 cm signal vanishes except in 
the high density neutral regions. Since the observations 
directly probe the neutral hydrogen density, one can use 
it to probe the detailed topology of the ionized regions in 



the pre-overlap phase. It is therefore essential to model 
the clustering of the sources accurately so as to predict 
the reionization topology. 



2. Models of source clustering 

There have been various approaches to account for the 
clustering of sources, most of them using the properties of 
the gaussian random field in terms of the extended Press- 
Schechter formalism. We have written equation in 
terms of globally-averaged quantities; now let us write it 
in a slightly different form where the averaging is done 
over a spherical region of radius R which has a density 
contrast 6 (linearly extrapolated to present epoch). Then 

Ti = ^ion t;^ {Qmi)6,RCa{T)nH(l + d) 

at at 

(51) 

where we have assumed that reionization is primarily 
driven by galaxies and have used equation to write 
riph in terms of /coii- The above equation can be solved 
for a given {5, R} if we know the form of (/co1i)i5,_r- It 
turns out that one has a simple generalization of equation 
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p5)) which encodes the clustering of sources at different 
scales and is given by js^ 



(/coll) 5,, 



erfc 



5jD{t) - 5 



V2[(T2(Af^i„)-a2(Af)] 



(52) 



where M = 47ri?'^p„/3 is the mass scale corresponding 
to R. The evolution of {Qmi)s,R can be used for study- 
ing the filling fraction of ionized regions within the IGM 
on various scales i? as a function of overdensity 5. In 
typical scenarios, this model predicts that reionization 
is driven by overlap of individual ionized regions around 
clustered sources residing in overdense regions of the Uni- 
verse. This leads to an "inside-out" scenario of reioniza- 
tion where, on average, high-density regions are ionized 
first. 

In a somewhat similar but slightly different approach, 
one can obtain the size distribution of the ionized regions 
at a given epoch. If we integrate equation ([5T|) upto a 
certain time i, we can write the filling factor as 



(Qhii) 



S.R 



coU)5,R 



1 + N,, 



(53) 



where (Qhii)^,^ A^rcc = J dt' {Qmi)s,R C a(r)nH(l + <5) 
is the number of recombinations within the region over 
the history of the IGM. Hence, the condition for the re- 
gion to be fully ionized (Qhtt)a « > 1 is g iven by a con- 
dition on the collapsed faction [82 



{fcon)sM>NrJ^{l + N,,,) 



(54) 



The condition for a region to be "self-ionized" can be con- 
verted into a condition in terms of the density contrast 
d. The problem then is very similar to the problem of 
collapse of haloes where it is studied whether the density 
averaged over a spherical volume exceeds a critical value 
( "barrier" ) [s^ . The only difference here is that the bar- 
rier is much complex than the collapse of halo problem. 
One can approximate the barrier with a linear one and 
then write a distribution of sizes similar to equation ([38|) . 
The results obtained from this approach is very similar to 
the one obtained from the previous one that reionization 
at large scales proceeds in a "inside-out" fashion. 

Both the approaches described above have also been 
extended to simulation boxes and used for making mock 
maps of the neutral hydrogen distribution which is ex- 
tremely useful for 21cm observations. These methods es- 
sentially give a semi-analytic (or semi-numeric) approach 
to deal with the radiative transfer problem and can be 
used for making maps using much less computing power. 



3. Recombination and self- shielding 

In the above approaches, recombinations in the IGM 
an be accounted for only by averaging over a spherical re- 
gion. In reality, even if a given spherical region contains 
enough photons to self-ionize, the high-density clumps 




FIG. 3: Ionization maps for a IGM with ionized hydrogen 
fraction of 75% for models with different assumptions regard- 
ing the spatial distribution of sinks and sources of ionizing 
radiation [83|. The black regions represent ionized ones. The 
left panel considers only homogenous recombination while 
the middle and right panels consider different levels of self- 
shielding. The thickness of the slice shown is l/i~^Mpc. 



within the region will remain neutral for a longer period 
because of their high recombination rate and thus alter 
the nature of the ionization field. A simple prescription 
to describe the presence of such neutral clumps by assum- 
ing that regions with overdensities above a critical value 
(A > Ac) remain neutral as we discussed earlier. Un- 
fortunately this is also not fully appropriate as many of 
the high-density regions are expected to harbour ionizing 
sources. Whether a region remains neutral will depend 
on two competing factors, the local density (which de- 
termines the recombination rate) and the proximity to 
ionizing sources (which determines the number of pho- 
tons available). It is thus important to include a realistic 
spatial distribution of recombinations into the formalisms 
for making ionization maps. 

A possible way of modelling the recombinations in 
high-density regions is to use a self-shielding criterion. In 
order to be ionized, a given point should satisfy the con- 
dition that it cannot remain self-shielded, i.e., it should 
not be a Lyman-limit system. In terms of the column 
density of the system, the condition can be written as 



A^Hi = nHi(x)L(x)CTHi < 1, 



(55) 



where nHi(x) is the number density of neutral hydro- 
gen at the given point and L{x.) is the size of the of 
the absorber. It turns out that the topology of reioniza- 
tion is very different if such self-shielding is taken into 
account particularly when reionization is extended and 
photon-starved as can be seen from Figure [31 Initially 
reionization proceeds inside-out with the high density re- 
gions hosting the sources of ionizing sources becoming 
ionized first. In the later stages, the high density regions 
which are far from sources remain neutral and reioniza- 
tion proceeds deep into the underdense regions before 
slowly evaporating denser regions. Such models can, in 
principle, be constrained by the first generation of 21 cm 
experiments. 
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IV. CURRENT STATUS AND FUTURE 



the IGM and hence extend the reionization process. 



In this section, let us review the current status of vari- 
ous approaches to studying reionization and their future 
prospects. 



A. Simulations 

Though the analytical studies mentioned above allow 
us to develop a good understanding of the different pro- 
cesses involved in reionization, they can take into account 
the physical processes only in some approximate sense. 
In fact, a detailed and complete description of reioniza- 
tion would require locating the ionizing sources, resolving 
the inhomogeneities in the IGM, following the scattering 
processes through detailed radiative transfer, and so on. 
Numerical simulations, in spite of their limitations, have 
been of immense importance in these areas. 

Since the ionizing photons during early stages of reion- 
ization mostly originate from smaller haloes which are 
far more numerous than the larger galaxies at high red- 
shifts. The need to resolve such small structures requires 
the simulation boxes to have high enough resolution. On 
the other hand, these ionizing sources were strongly clus- 
tered at high redshifts and, as a consequence, the ionized 
regions they created are expected to overlap and grow to 
very large sizes, reaching upto tens of Mpc [88,, ^90]. 
As already discussed, the many orders of magnitude dif- 
ference between these length scales demand extremely 
high computing power from any simulations designed to 
study early structure formation from the point of view of 
reionization. 

To simulate reionization, one usually runs a N-body 
simulation (either dark matter only or including baryons) 
to generate the large-scale density field, identifies haloes 
within the density field and assign ionizing photons to the 
haloes using a assumption like equation (j48p . It turns out 
that the most difficult step is to solve the radiative trans- 
fer equation and study the growth of ionized regions. In 
principle, one could solve equation ^ directly for the in- 
tensity at every point in the seven-dimensional {t, x, n, i/) 
space, given the absorption coefficient and the emissivity. 
However, the high dimensionality of the problem makes 
the solution of the complete radiative transfer equation 
well beyond our capabilities, particularly since we do not 
have any obvious symmetries in the problem and often 
need high spatial and angular resolution in cosmologi- 
cal simulations. Hence, the approach to the problem has 
been to use different numerical schemes and approxima- 
tions, like ray-tracing [T^. [oi. lO' 



931194 19a 1961. 1971, |98t, 

10ll| , local depth approx- 
At present, most of the 



Monte C arlo methods [99', lOCf, 
imation jl02| and others [103]. 
simulations do not have enough resolution to reliably 
identify the low mass ^ IO^Mq sources which were prob- 
ably responsible for early stages of reionization. Also, 
there are difficulties in resolving the small scale struc- 
tures which contribute significantly to the dumpiness in 



B. Various observational probes 

Finally, we review certain observations which shape 
our understanding of reionization. 



1. Absorption spectra of high redshift sources 

We have already discussed that the primary evidence 
for reionization comes from absorption spectra of quasars 
(Lya forest) at z < 6. We have also discussed that the 
effective optical depth of Lya photons becomes signifi- 
cantly large at z ^ 6 implying regions with high trans- 
mission in the Lya forest becoming rare at high redshifts 
[35I HI]. Therefore the standard methods of analyzing 
the Lya forest (like the probability distribution function 
and power spectrum) are not very effective. Amongst 
alternate r nethqds, one can use the the distribution of 
dark gaps [l04l llOSi ] which are defined as contiguous re- 
gions of the spectru m ha ving an optical depth above a 
threshold value [s^, llOSi ]. It has been found that the 
current observations constr ain t he neutral hydrogen frac- 
tion xm < 0.36 at 2: = 6.3 [l06| ]. It is expected that the 
SDSS and Palomar-Quest survey Il07| would detect ~ 30 
quasars at these redshifts within the next few years and 
hence we expect robust conclusions from such studies in 
very near future. 

Like quasars, one can also use absorption spectra of 
other high redshift energetic sources like gamma ray 
bursts (GRBs) and supernovae. In fact, analyses using 
the damping wing effects of the Voigt profile have been 
already performed on the GRB detected at a redshift 
z = 6.3, and the wing shape is well-fit by a neutral frac- 
tion xm < 0.17 jl08j ]. The dark gap width distrib ution 
gives a similar constraint xhi — (6.4 ± 0.3) x 10~^ [l09| ]. 
In order to obtain more stringent limits on reionization, 
it is important to increase the sample size of z > 6 GRBs. 

In addition to hydrogen reionization, the Lya forest in 
the quasar absorption lines at 2; w 3 can also be used 
for studying reionization of singly ionized helium to dou- 
bly ionized state (the reionization of neutral helium to 
singly ionized state follows hydrogen for almost all types 
of sources). The helium reionization coincides with the 
rise in quasar population at z '--^ 3 and it effects the ther- 
mal history of the IGM at these redshifts. However, there 
are various aspects of the observation that are not well 
understood and requires much deta iled modelling of he- 
lium reionization [UOi iuA EH • 



2. CMBR observations 

As we have discussed already, the first evidence for an 
early reionization epoch came from the CMBR polariza- 
tion data. This data is going to be much more precise in 
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future with experiments like PLANCK, and is expected 
to improve the constraints on t^i. With improved statisti- 
cal errors, it might be possible to distinguish between dif- 
ferent evolutions of the ionized fraction, particularly with 
E-mode polarization a uto-correla tion, as is found from 
theoretical calculations [ll4l Ill5j | . An alternative option 
to probe reionization through CMBR is through the small 
scale observations of temperature anisotropics. It has 
been well known that the scattering of the CMBR pho- 
tons by the bulk motion of the electrons in clusters gives 
rise to a signal at large multipoles £ ^ 1000, known as the 
kinetic Sunyaev Zeldovich (SZ) effect. Such a signal can 
also originate from the fluctuations in the distribution of 
free electrons arising from cosmic reionization. It turns 
out that for reionization, the signal is dominated by the 
patchiness in the ng-distribution. Now, in most scenarios 
of reionization, it is expected that the distribution of neu- 
tral hydrogen would be quite patchy in the pre-overlap 
era, with the ionized hydrogen mostly contained within 
isolated bubbles. The amplitude of this signal is signif- 
icant around £ ~ 1000 and is usually comparable to or 
greater than the signal arising from standard kinetic SZ 
effect. Theoretical estimates of the signal have been per- 
formed for various reionization scenarios, and it has been 
predicted that the exper iment can be used for constrain- 
ing reionization history jll6l . Also, it is possible to 
have an idea about the nature of reionization sources, as 
the signal from UV sources. X-ray sources and decaying 
particles are quite different. With multi-frequency exper- 
iments like Atacama Cosmolog y Tel escope (ACT) |l54l | 
and South Pole Telescope (SPT) ji55| coming up m near 
future, this promises to put strong constraints on the 
reionization scenarios. 



Lya emitters at high redshifts is complicated by various 
factors like the velocity of the sources with respect to the 
surrounding IGM, the density distribution and the size of 
ionized regions around the sources and the clustering of 
sources. It is thus extremely important to have detailed 
models of Lya emitting galaxies in order to use them for 
constraining reionization. 



4-. Sources of reionization 

As we discussed earlier, a major challenge in our under- 
standing of reionization depends on our knowledge of the 
sources, particularly at high redshifts. As we understand 
at present, neither the bright z > 6 quasars discovered 
by the SDSS gro up [l24| nor the faint ones detected in X- 
ray observations |125l | produce enough photons to reion- 
ize the I GM . The disc overy of star-forming galaxies at 
2 > 6.5 [HE EH Il28j | has resulted in speculation that 
early galaxies produce bulk of the ionizing photons for 
reionization. Unfortunately, there are significant uncer- 
tainties in constraining the amount of ionizing radiation 
at these redshifts because the bulk of ionizing photons 
could be produced by faint sources which are beyond the 
present sensitivities. In fact, some models have predicted 
that the z > 7 sources identified in these surveys are rel- 
atively massive (M « IO^A/q) and rare objects which 
are only marginally (f* 1% ) contributing to the reion- 
ization photon budget [l29| |. A much better prospect of 
detecting these sources would be through the Ultra-Deep 
Imaging Survey using the future telescope JWST. 



3. Lya emitters 

In recent years, a number of groups have studied star- 
forming galaxies at z ~ 6 — 7, and measurements of the 
Lya emission line luminosity function evolution provide 
another useful observational constraint [nlUTg^. While 
the quasar absorption spectra probe the neutral hydro- 
gen fraction regime xm < 0.01, this method is sensitive 
to the range xm ~ 0.1 — 1.0. Lya emission from galaxies 
is expected to be suppressed at redshifts beyond reioniza- 
tion because of the absorption due to neutral hydrogen, 
which clearly affects the evolution of the l uminosity func- 
tion of such Lya emitters at high redshifts |ll8l . ll20lll2l| . 
Thus a comparison of the luminosity functions at differ- 
ent redshifts could be used for constraining the reioniza- 
tion. Through a simple analysis, it was found that the 
luminosity functions at z — 5.7 and z = 6.5 are statis- 
tically consistent with one another thus implying that 
reionization was largely complete at z « 6.5. More so- 
phisticated calculations on t he evolution of the luminos- 
ity function of Lya emitters |ll8l . Il20l . Il22j | suggest that 
the neutral fracti on o f hydrogen at z = 6.5 should be less 
than 50 per cent [l23| | . Unfortunately, the analysis of the 



5. 21cm experiments 

We have already discussed the basic theory behind 
detecting the fluctuations in the neutral hydrogen den- 
sity through the 21 cm emission experiments. There are 
essentially two complementary approaches to studying 
reionization using 21 cm signal. The first one is through 
global statistical properties of the neutr a l hydrogen sig- 
nal, like the power spectrum jl30l . Il3ll Il32l Il33l . Il34 
Il35l |. The second one is to directly detect the ionized 
bubbles around sources, either through blind surveys or 
via targetted observations [HI, [Hi, [iM [HI, [TH, [mO] . 

The major difficulty in obtaining the cosmological sig- 
nal from these experiments is that it is expected to be 
only a small contribution buried deep in the emission 
from other ast roph ysical sour ces (foregr ounds ] and in the 
system noise [mI [TH [TH, [HI [mI [HI tHl- It is 
thus a big challenge to detect the signal which is of cos- 
mological importance from the other contributions that 
are orders of magnitude larger. Once such challenges are 
dealt with, this probe will be the strongest probe for not 
only reionization, but of the matter distribution at very 
small scales during the dark ages. 
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V. CONCLUDING REMARKS 

We have discussed the analytical approaches to model 
different aspects of rcionization which will help in un- 
derstanding the most relevant physical processes. In an 
explicit example, we have shown how to apply this for- 
malism for constraining the reionization history using a 
variety of observational data. These constraints imply 
that reionization is an extended process over a redshift 
range 15 > 2; > 6. It is most likely driven by the first 
sources which form in small mass haloes. However, there 



are still uncertainties about the exact nature of these 
sources and the detailed topology of ionized regions. Such 
details are going to be addressed in near future as new ob- 
servations, both space-borne and ground-based, are likely 
to settle these long-standing questions. From the theo- 
retical point of view, it is thereby important to develop 
detailed analytical and numerical models to extract the 
maximum information about the physical processes rele- 
vant for rcionization out of the expected large and com- 
plex data sets. 
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